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“Algebraic  and  Computational  Aspects  of  Network  Reliability  Problems” 


SUMMARY 

This  research  has  advanced  both  theoretical  and  computational  aspects  of 
evaluating  the  reliability  of  a  complex  system  in  terms  of  its  structure  and  the 
reliability  of  its  components.  This  type  of  problem  arises  in  particular  in  the  design 
and  evaluation  of  telecommunication  and  distribution  systems,  which  are  commonly 
modelled  as  networks.  The  present  research  employs  an  algebraic  approach  for 
studying  the  reliability  of  such  network  systems.  This  approach  has  not  only  unified 
certain  theoretical  aspects  of  network  reliability  problems  but  has  also  suggested  a 
number  of  new  algorithms  for  calculating  various  reliability  measures.  Based  on  this 
approach,  both  exact  and  approximate  computational  schemes  have  been  devel¬ 
oped,  together  with  supporting  data  structures  for  implementing  the  necessary 
computations  in  an  efficient  manner.  Approximation  schemes,  also  based  on  an 
underlying  algebraic  structure,  have  also  been  developed  for  evaluating  more 
general  measures  of  system  performance,  such  as  average  delay  or  throughput. 

RESEARCH  OBJECTIVES 

The  objective  of  this  research  project  is  to  study  network  reliability  problems 
by  identifying  and  exploiting  certain  underlying  algebraic  structures.  These 
structures  enable  the  exact  calculation  of  reliability  using  reasonably  efficient 
(“pseudopolynomial”)  algorithms  applicable  to  certain  special  classes  of  networks. 
The  underlying  structure  governing  this  situation  is  that  of  a  semilattice,  defined  on 
certain  combinatorial  objects  (paths  and  cutsets)  and  for  which  a  general  recursive 
algorithm  can  be  used  for  calculating  system  reliability.  In  addition,  approximation  of 
the  reliability,  or  more  general  measures  of  system  performance,  derives  from 
studying  a  different  type  of  algebraic  structure  (a  distributive  lattice)  in  more  general 
network  systems. 

ACCOMPLISHMENTS 

The  calculation  of  network  reliability  in  general  is  known  to  be  mathematically 
difficult  (NP-hard).  Thus,  a  more  promising  line  of  attack  is  to  investigate  the 
calculation  of  reliability  when  the  underlying  network  has  certain  special  structure. 
One  commonly  occurring  structure  is  that  of  a  planar  network:  a  network  that  can  be 
represented  in  the  plane  without  crossing  edges.  An  algebraic  structure  (a  lattice) 
governing  planar  networks  has  been  identified  and  studied.  This  algebraic  structure 
enables  the  “pseudopolynomial”  calculation  of  network  reliability;  that  is,  the  compu- 
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tational  effort  is  polynomially  bounded  in  terms  of  the  number  of  paths  (or  cutsets) 
joining  two  specified  nodes  in  the  network.  To  streamline  these  calculations,  efficient 
techniques  for  generating  paths  and  cutsets  in  planar  networks  have  also  been 
developed.  The  resulting  algorithm  has  been  developed  and  applied  to  a  number  of 
difficult  benchmark  problems  in  the  literature,  providing  in  many  cases  the  first 
complete  reliability  analysis  of  such  problems.  In  one  particular  case,  the  compu¬ 
tation  time  was  reduced  from  the  order  of  hours  to  the  order  of  minutes.  Another 
feature  of  the  algebraic  approach  has  been  that  it  enables  the  calculation  of  the 
entire  reliability  polynomial  in  functional  form,  rather  than  simply  the  numerical  value 
at  a  particular  setting  for  the  input  edge  reliabilities.  This  research,  conducted  jointly 
by  D.  Shier,  J.  P.  Jarvis,  and  D.  Whited  (Ph.D.  student),  has  resulted  in  a  number  of 
papers  [3.  6, 1 1 , 1 2]  as  well  as  the  Ph.D.  dissertation  [1 0].  Recent  collaboration  with 
M.  Ball  and  J.  S.  Provan  has  extended  these  ideas  to  the  pseudopolynomial  calcu¬ 
lation  of  reliability  for  other  types  of  coherent  systems;  generalizations  of  consecutive 
k-out-of-n  systems  and  reliability  covering  problems. 


An  algebraic  approach  has  also  led  to  several  new  methods  for  approximating 
the  s-t  reliability  of  a  network  (the  probability  that  an  operative  path  exists  between 
two  given  nodes  s  and  t  in  the  network).  The  method  developed  is  an  iterative  one, 
which  produces  both  a  sequence  of  lower  bounds  and  a  sequence  of  upper  bounds 
on  the  exact  s-t  reliability  value.  It  has  been  proved  that  both  the  lower  and  the  upper 
bounds  will  converge  to  the  correct  answer  after  a  finite  number  of  iterations.  More¬ 
over,  the  lower  bounds  produce  a  nondecreasing  sequence  of  polynomials  and  the 
upper  bounds  produce  a  nonincreasing  sequence  of  polynomials,  converging  to  the 
true  reliability  polynomial.  Thus  at  each  step  one  has  an  interval  guaranteed  to 
contain  the  correct  value;  it  is  then  possible  to  terminate  the  iterative  process  once 
this  interval  is  sufficiently  small.  This  aspect  is  particularly  important  since  the  exact 
calculation  of  reliability  can  be  quite  time-consuming,  even  for  moderately  small 
networks.  The  algebraic  approach  developed  here  produces  a  general  polynomial 
expression  for  the  s-t  reliability  in  terms  of  the  input  parameters  rather  than  simply  a 
single  numerical  answer.  This  symbolic  form  is  especially  useful  for  assessing  the 
“importance”  of  individual  links  in  the  original  network,  as  well  as  facilitating  a 
general  sensitivity  analysis.  These  results,  obtained  by  D.  Shier  and  D.  Whited,  are 
discussed  in  the  papers  [7,  8]. 

Another  polynomial,  related  to  the  reliability  polynomial,  is  the  chromatic  r 
polynomial  of  an  undirected  network  whose  coefficients  reveal  a  good  deal  of  infor-  ^ 
mation  about  the  combinatorial  structure  of  the  network.  In  particular,  one  of  these  Q. 

coefficients  gives  the  “all  terminal  domination”  of  the  network,  an  invariant  that  has  □ 

been  recently  studied  in  some  detail  in  relation  to  the  “all  terminal  reliability”  of  a  ■  ■ 
network  (the  probability  that  all  nodes  of  the  network  can  communicate).  Our 
research  has  provided  more  efficient  methods  for  calculating  the  chromatic  poly-  ■  ■" 
nomial,  based  on  decomposition  of  the  original  network  into  trees  [1]  and  into  - 
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chordal  subgraphs  [4].  This  algorithmic  research,  conducted  jointly  by  D.  R.  Shier 
and  graduate  students  S.  Frank  and  N.  Chandrasekharan,  is  reported  in  [1 , 4]. 

The  algebraic  focus  of  this  research  has  also  been  applied  to  the  calculation 
of  general  performance  measures  (including  reliability,  delay,  throughput)  defined  on 
a  network  with  stochastically  failing  components.  Again,  since  the  exact  calculation 
of  such  measures  is  known  to  be  mathematically  intractable,  we  have  instead  devel¬ 
oped  approaches  for  approximating  the  performance  of  the  system,  based  on  gener¬ 
ating  a  relatively  small  subset  of  the  “states”  of  the  system  that  covers  in  probability  a 
large  proportion  of  the  total  probability.  It  has  been  found  that  an  elegant  algebraic 
structure  (a  distributive  lattice)  governs  the  case  in  which  each  component  assumes 
one  of  two  modes:  “working"  or  “failed.”  This  lattice  turns  out  to  have  a  number  of 
significant  properties;  it  is  rank  unimodal,  rank  symmetric,  and  satisfies  the  Sperner 
property.  Study  of  the  theoretical  properties  of  this  lattice  has  yielded  a  reasonably 
effective  algorithm  for  generating  in  order  the  most  probable  states  of  the  system, 
without  the  need  to  examine  the  entire  state  space.  In  one  fairly  typical  example 
having  25  components,  only  0.08%  of  the  states  of  the  system  needed  to  be  gen¬ 
erated  to  cover  95%  of  the  total  probability.  Rather  interestingly,  the  computational 
complexity  of  the  new  algorithm  turns  out  to  be  related  to  a  certain  combinatorial 
invariant  of  this  lattice.  A  number  of  papers,  written  jointly  by  D.  Shier,  R.  Jamison, 
and  E.  Valvo  (graduate  student),  have  been  devoted  to  discussing  this  approach  [2, 
5,  9].  More  recent  work,  involving  J.  P.  Jarvis,  D.  Shier,  R.  Lakin  (graduate  student), 
and  E.  Bibelnieks  (graduate  student),  has  considered  the  more  general  case  in 
which  each  component  of  the  system  can  assume  a  finite  number  of  states.  Such  a 
model  is  appropriate  when  a  system  degrades  slowly  over  time,  rather  than  simply 
fails.  It  appears  that  an  algebraic  approach  will  also  be  successful  in  this  situation. 
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